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Multi-pion correlations and wavepacket size effects on the 
pion multiplicity distribution, pion momentum distribution 
and two-pion interferometry are studied. It is shown that 
multi-pion Bose-Einstein correlations and the wavepacket size 
cause an abundance of pions at low momentum, increase the 
mean pion multiplicity and decrease both the apparent radius 
of the source and the coherent source parameter derived from 
two-pion interferometry. 

PACS number(s): 25.75 Gz,11.38 Mh, 11.30 Rd. 



I. INTRODUCTION 

The main goal of colliding heavy nuclei at high energy 
is to produce quark-gluon plasma (QGP). One of the key 
quantities related to QGP is the energy density. To deter- 
mine the energy density, we need to know the space-time 
information of the source which is the main objective 
of Hanbury-Brown-Twiss correlations [Q,^]. Two-pion 
Bose-Einstein(BE) correlations are widely used in high 
energy heavy-ion collisions to provide information on the 
space-time structure, degree of coherence and dynamics 
of the region where the pions were produced . 

In principle, the extension of two-pion interferometry 
to multi-pion interferometry is straightforward [p|-p3[. 
Experimentally, ultrarelativistic hadronic and nuclear 
collisions provide the environment for creating dozens, 
and in some cases hundreds, of pions [p4j-p7|. Therefore, 
it is necessary to consider multi-pion correlations in these 
processes. The bosonic nature of the pion should affect 
the single and two-pion spectra and distort the two-pion 
correlation function. There is a kind of coherence length 
which corresponds to the wave packet length of the emit- 
ter, which cause pions to concentrate at low momenta 
pq fjof . Thus, it is very interesting to analyse the effects 
of the multi-pion correlations and wavepacket length on 
these variables. 

The arrangement of this paper is as follows: The 
derivation of the multi-pion correlation function for a to- 
tally chaotic source is given in section 2. In section 3, 
the effects of multi-pion correlations on the multiplicity 
distribution are discussed. In section 4, the effects of 
multi-pion correlations on two-pion correlation function 
are analysed for n-pion events and for events with all pos- 
sible multiplicities. In section 5, our results for a source 
distribution are given. Conclusions are given in section 
6. 



II. PION CORRELATION FUNCTION FOR A 
CHAOTIC SOURCE 

The general definition of the pure n pion correlation 
function C n (pi, • • -p n ) is 



C n (pi, ■ ■ -Pn) 



PniPl 



■Pn 



(i) 



where P n (pi , • • -p n ) is the probability of observing n-pions 
with momenta {p{\ all in the same n-pion event. 

A state created by a classical pion source is described 

by imn 



\4> >= exp(i J dp J d A xj{x) exp(ipx)c + (p)\0 >, (2) 

where c + (p) is the pion-creation operator. j(x) is the 
pion current, which can be expressed as 



j(x) = / d A x' d A pj(x' 1 p)v(x') exp(— ip(x — x')). (3) 



Here j(x',p) is the probability amplitude of finding a 
pion with momentum p , emitted by the emitter at x' . 
v(x') = exp(i<fi(x')) is a random phase factor which has 
been extracted from j(x',p). All emitters are uncorre- 
cted in coordinate space when assuming 

{v*{x')v{x)} = 5 i {x'-x). (4) 

Here {• • •} means phase average which is defined as j|] 
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(5) 



Eq.(4) is the ideal case. In more realistic cases, each 
chaotic emitter has a small coherent wave packet length, 
and Eq. (4) can be replaced by 



(x 3 - x' 3 ) 2 (x - x' ) 



' \2 



s 2 



s 2 



(6) 



Here 5 is a parameter which determines the coherent 
length (time) scale of the emitter^]. For simplicity, the 



1 S is smaller than the total size (R) of the source which 
consists of all emitters 
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same coherent scale is taken for both space and time at Here Y — Xl + X2 and y = x% — x 2 are four dimensional 

the moment. The above formula shows that two emit- coordinates, while k = Pl + P2 and q = p\ ~p 2 are the cor- 
ters within the range of 5 can be seen as one emitter, responding four dimensional momenta. g w (Y,k) is the 
while two-emitters out of this range are incoherent. From Wigner function , which can be understood as the prob- 
Eq.(5) we have ability of finding a pion at Y with momentum k. It is 

defined as 



j'=i 



(7) 



i=l 



The coherent state can be expanded in Fock-Space as 
(i J j{x)e lpx c + {p)dpd 4 x) 



n=0 



|0 >=£>>, (8) 

n=0 



with 



^ >= (j J j(x)e i P x c+(p)dpdx) n [Q > 



(9) 



Here \n > is an rt-pion state. According to the definition 
of Eq.(l), the pure n-pion correlation function can be 
re-expressed as 

C n {p\, ■ ■ ■p n )= 

{< n\c + (pi) ■ ■ ■ c + {p n )c(p n ) ■ ■ ■ c(pi)\n >} 

nr=i{<iic + (Pi)c(Pi)ii>} 

From above definition, the pure two-pion correlation 
function reads (see Appendix) 



g w (Y, k) = {J d 4 yf(Y + y/2)j(Y - y/2) exp(-ik ■ y)}. 

(14) 

Inserting Eq.(3) into the above equation we have 

d^x 1 j*{x' iPx^pxe^^+y^-^v*^ 1 ) (15) 
d i x"j(x",p 2 )d i p 2 e- t P^ Y -y/ 2 - x '">v{x")}, 



so the pure single pion inclusive distribution P\{p) can 
be expressed as: 



Px(k) = Jg w (Y,k)d 4 Y 



(10) 



C 2 (pi,P 2 ) 



{< 2|c + (pi)c+(p 2 )c(^)c(pi)|2 >} 
{< l|c+(pi)c(pi)|l >}{< l\c+(p 2 )c(p 2 \l >} 

= 1 + 

{< ljc+(gi)c(g^l >}{< l|c+(pb)c(pi)|l >} 
{< l|c+(ft)c(ft)|l >}{< l\c+(p 2 )c(p 2 )\l >}■ 

Using the relationship 



(11) 

with fco 
reads 



{jd 4 x'd 4 x"d 4 p 1 d 4 p 2 r(x', Pl )j(x",p 2 ) 

u*{x')v{x")8 4 {k - p ^)S 4 ( Pl - pa) 

e ipi(Y-x') e -ip 2 (Y-x")^ 

{Jd 4 x'd 4 x"j*(x',k)j(x",k) 
v*(x')v(x")e- ik< ~- x '- x "^} 



(16) 



k 2 + m%. The two-pion correlation function 



c(p)\n >= i I d 4 xj(x) exp(ipx)\n — 1 > 
we have 



(12) 



C 2 {pi,p 2 ) = 1 + 



J d 4 xd 4 x'g w (x, k)g w (x', k) expjiq ■ (x - x')) 
J d 4 xd i x'g w (x,p 1 )g w (x' ,p 2 ) 

Similarly, the pure n-pion correlation function can be 
expressed as 

C„(pi,p 2 , ' ' -,Pn) = £xi,a(l)X2, < r(2)-Xn, -(n) (18) 

with 

Xij=x(Pi, Pi )= J 9wil 2 } (19) 

y J d 4 xd i yg w (x, p t )g w (y , p ■) 

2 I would like to express my thanks to the referee and Dr. Here a(i) denotes the ith element of a permutation of 
Hans Feldmeier who make this definition more clear. the sequence 1, 2, 3, • • •, n, and the sum over a denotes 

the sum over all n\ permutations of this sequence. 



{< l|c+(ft) C (fc)|l >} = {J d 4 x 1 d 4 x 2 f(x 1 )j(x 2 ) 

exp(-i(pixi -p 2 x 2 ))} 
= {J d 4 yd 4 Y.f(Y + y/2) 
j{Y - y/2) exp(-iky - iqY)} 
d 4 Yg w (Y, k) exp(-iq ■ Y). (13) 
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In the following, we shall often use the expression of the 
pure n-pion momentum probability distribution P n (pi, • • 
■p n ) given by 

Pn(fl,- ■ ;Pn) = ^Pl,a(l)P2,a(2)---Pn,a(n) (20) 



with 



1 [ ,4 , (Pi+Pi)- 



Pid = PiPi'Pj) = — / d x 9w(x 



??.() 



J{pi-Pj)x 



(21) 



Here no is the mean multiplicity without BE correla- 
tions, defined by 



n o = / g w {x,k)d A xdk, 



(22) 



Once the source distribution g w (x,k) is known, the 
pure multi-pion correlation function can be calculated 
from Eq.(18). 



III. EFFECTS OF MULTI-PION 
CORRELATIONS ON PION MULTIPLICITY 
DISTRIBUTION 

From the above definition of |n > we have 
w(n) = {< n\n >} = ^ J P n (p u ■ ■ -,p n ) f[ dpi (23) 

i— 1 

and 

{<</#>} = ]T{< fc|£; >} 
k 

n k r 

= ^2 "w / ' ' ' d Pk p k(pi, ■ ■ -,Pk) 

k ' ^ 

= £"(fc). ( M ) 
fc 

Then the pion-multiplicity distribution becomes 
_ {< n|n >} _ iojn) 

P{n) ~ {< 0|0 ( } 

For particles without BE correlations, we have 

pi,j = $i,jPi,i (26) 

and 



J Pn(Pl, ■ ■ -PrOYidPi = !' 



therefore 



P(n) = e 



„-no 



(27) 



(28) 



That is P(n) reduces to the usual Poisson distribution. 
If we take into account multi-pion Bose-Einstein correla- 
tions, the pion multiplicity distribution is not a Poisson 
distribution anymore (see Fig.l). From the above equa- 
tion, we can also obtain the mean multiplicity < M > 



< M >= V n ■ P{n) 



(29) 



IV. EFFECTS OF MULTI-PION CORRELATIONS 
ON TWO-PION CORRELATION FUNCTION 

A. Two-pion correlation function for n-pion events 

For n(n > 2)n events, the two-pion correlation func- 
tion can be defined as [] 



C 2 n (f>i,p 2 ) 



i?(P!,P2) 
P{ 1 {PI)P{ 1 (P2Y 



(30) 



where P"(px, ■ ■ -pi) is the modified i-pion inclusive distri- 
bution in n pion events. The definition of P™(pi, • • -,Pi) 

isFI 



{< n|c+(pi) ■ ■ • c+(pi)c(pi) ■ • ■ c(pi)\n >} 



{< n\n >}n • ■ ■ (n — i + 1) 



(31) 



For an event with multiplicity n, the modified two- 
pion inclusive and single-pion inclusive distributions can 
be expressed as 



and 



pn (r j IlTU d P^ P n(Pl--Pn) 
J 1 li=l dPiPn{Pl--Pn) 

JTVi=2 d Pi P n(Pl-~-Pn) 



/nr=i dp l p n (p 1 ....p n ) 



(32) 



(33) 



respectively. 

As n increases, the computation of the above intcrgrals 
becomes more and more complex. Now we define the 
function 



G i {p,q)=n l J p(p,pi)dpip(piip 2 )dp 2 

■ ■ ■/o(Pi-2,Pi-i)rfp'i-ipfe-i,?)- (34) 

From the expression of P n (pi> ■ ■ ■,Pn) (Eq. (20)), the 
two-pion inclusive distribution can be expressed as 



This definition takes into account the higher order Bose- 
Einstein correlations. 
4 It is easily checked that J dpi ■ ■ ■ dpiP"(pi, •••,$) = 1. 
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1 1 n 

P?(p,q) = - f TT^E 

n(n — 1) oj(n) 



E G m (p,p)G t - m (q,q) + 

m— 1 

G m (p, 9) ■ Gi_ m (9,p)]tj(n - i) (35) 



with 



win 



n /• n 

^ / II d PkPn(Pl, ■ ■ -,Pn 
J k=l 



The single-pion distribution is given as 
1 1 - 

fr(p1 = ~nEGifepM»-i). 

n uj(n) * — ' 
v ; i=i 



(36) 



(37) 



Now the main problem is to extract the expression of 
to(n). From Eq.(37) we have 



{<4>\c+{p)c+{q)c{q)c{p)\4>>} 

{< 4>\4> >} < M(M - 1) > 
{< <j>\<j> >} < M(M- 1) > 

1 y 

1) > z - 



{< ^\<j> >} < M (M 

n i—X 

G m{p,p)Gi- m (q,q) + 

i=2 m=l 

G rn (p, q) ■ G l - m (q,p))uj(n - i)} 

Ei,j gj(P.P)Gj(g, g) + Gi(p, g) ■ gj(g,p) 

< M(M- 1) > ' (j 

The two-pion correlation function for all multiplicity 
distribution is 



Pt(p,q) 
Pt{q)Pf{p)' 



(42) 



1 - 



with 



C(i) = / dpG t (p,p). 



(38) 



(39) 



Using the above method, the two-pion and single 
pion inclusive distributions can be calculated for n-pion 
events. Gi(p, q) and C(i) can be calculated using Monte- 
Carlo or analytical integration. 



B. Two-pion correlation function for all events 

For a state \<j> > which contains all possible multiplic- 
ities, the modified single pion spectrum distribution can 
be expressed as^] 

p* n = {<<p\c + {p) C (m>} 

1 [P) {< </>\</> >} < M > 
_ E n Pr(p)-n-Qj(n) 
{< 4>\(f> >} < M > 

<M> 



(40) 



{< (f)\(f> >}<M > 



where Eq.(24) is applied. The two-pion inclusive distri- 
bution can be expressed as 



5 In calculating Pf we mix all the multiplicities while in 
calculating P™ we fix the pion multiplicity as n. 



V. MULTI-PION CORRELATIONS FOR A 
CHAOTIC SOURCE 

In this section, we will give an example to investigate 
the wave packet length and the multi-pion correlations 
effects on pion multiplicity, single pion and two-pion dis- 
tributions. 

We assume that the chaotic emitter amplitude distri- 
bution is 

222 
j(x,k) = % M)exp( 1 I -)S(x ) 



exp(- 



2A 2 



(43) 



Where i?o and A are the parameters which represents 
the radius of the chaotic source and the momentum range 
of pions respectively. Bringing Eq.(6) and Eq.(43) into 
Eqs. (15-17), we can easily get the source distribution 
function g w (x, K) 



1 3 x 2 
g w (x,K) = n ■ {—^r) 2 exp(-— T )S(x ) 

(^) l exp{-^ } , 



(44) 



the normalized pure pion single particle spectrum distri- 
bution 



^(p) = (i)*e*p{-4}, 
and the pure two-pion correlation function 



G 2 (pi,P2) = 1 + exp{- 



g 2 4i? 4 
2 S 2 + 4R 2 1 



(45) 



(46) 
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with 



R% 



111 

j 

A 2 ' a2 a2 



R 2 s 2 



a 2 s 2 + ml 



(47) 



In the above formula, the definition of q is q — pi — P2 ■ 
Then the derived radius through pure two-pion interfer- 
ometry is 



R 2 



o2 4i? 2 

°6 2 + 4R 2 



(48) 



It is clear that as the coherence length S increases, 
the mean momentum of pions (Ao) becomes smaller 
(Eq.(47)). It can be seen clearly from Eq.(48) that as the 
coherence length 5 increases, the apparent radius derived 
from pure two-pion interferometry also becomes smaller. 

g w (x, £y£) can be expressed as 



p + q 



) = n Q 



(7^)3/3 



(^A 2 )3/2 



(49) 



then 



p{p,q) = — I 9w{x, 
no J 



P + 1 j c i(p-Q)o 



dx 



(2ttA 2 )3/2 



2 n2 {p + q)' 2 

e o . 



(50) 



Defining 

Gn(p,q)=n l p(p,pi) Y[ dpip(pi,p i+ i) 

J i=l,n-2 
dp n ^ip(p n ^ 1 ,q) 

and using Eq.(50), we easily obtain 

-a n (p 2 +q 2 )+g n p-q 



G n {p,q) = riQ ■ a n e" 



where 



9l 



a n +i = ai — — — , b n = a n + a\, g n+ i 

40 n 



ffnffl 

2b n ' 



and 



with 



1 , 1 



(51) 



(52) 



(53) 



(54) 



R% 1 

Ol = -T- + 



4 ' 4A 2 ' 



9i 



R% 1 



2 2A 2 



, ai 



^A 2 )3/ 2 ' 



(55) 



Now we can calculate the pion multiplicity, single par- 
ticle and two-particle pion distributions according to the 



formula given in the preceding sections. The pion multi- 
plicity distribution is shown in Fig.l. It is clear that the 
probability with high pion multiplicity is greater than 
that for a conventional Poisson distribution. This is con- 
sistent with the nature of bosons, that is, the emission 
of bosons encourages the emission of more bosons. We 
also find that as the radius and temperature decrease, the 
effect of BE correlations on the multiplicity distribution 
becomes larger. This can be explained by the fact that as 
the radius and the temperature decrease the quantum ef- 
fect become larger. It is also clear that as the wavepacket 
8 increases, the effect of BE correlations on the multiplic- 
ity distribution becomes larger. It can be seen clearly 
from Eq.(47) that as 8 increases, A (the effective tem- 
perature) becomes smaller therefore the quantum effect 
becomes larger. 
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FIG. 1. Pion multiplicity distribution with and without 
(solid cure) Bose-Einstein correlations. The input value 
of A and no are 0.25GeV and 13. The dashed, dotted 
and dot-dashed line corresponds to Ro = 3/m ,5 = 0/m; 
Ro = 3/m ,5 — 0.5/m and Ro = 2.8/m ,5 = 0/m respec- 
tively. 

In Fig. 2, the mean multiplicity < M > is plotted 
against the parameter no, the mean multiplicity with- 
out BE correlations. When n is not very large, < M > 
increases with increasing no. However, at a critical point 
n c , the mean multiplicity < M > becomes to be infin- 
ity. From Eqs. (39-40), the mean multiplicity can be 
expressed as: 



< M >=£c(i) = y,« 



2a,i - g l 



x3/2 



(56) 



It is clear that the above equation will remain conver- 
gence provided that 



lim 

n— >oc 



C(n+1) 
C(n) 



_lim n • (^ ) 



(f) 3/2 ( 

On 

From Eqs. (53-55), we have 



1 

■ a 2 
2a n +i 



9n+l 



2a r , 



f /2 < 1. (57) 



5 



lim a n = 



Rg 
2A ' 
Rg I 



lim b n = ( h , . 

n^oo 2 2Aq 



lim g n = 0. 



(58) 
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FIG. 2. TVie mean multiplicity < M > t»s. no. T/ie 
soKrf ime corresponds to the case without BE correlation 
(no —< M >). The dashed and dotted lines corresponds to 
8 = Ofm and 8 = 0.5 fm respectively. The input value of Ro 
and A are 5/m and 0.36GeV . 

Eq.(56) remains convergence if 
n < l(l + R G A ) 3 




4i?2A 2 



4+(i?2A 2 + l)^ 2 /i?2- 



(59) 



We denote this critical pion multiplicity as n c . It is very 
interesting to notice that the critical multiplicity n c de- 
pends only on RqA.q and 5/Rq. The critical multiplicity 
n c vs. RoAq for different values of 5/Rq are shown in 
Fig. 3. It is clear that as 5/Rq increases the correspond- 
ing critical multiplicity becomes smaller. 
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FIG. 3. The critical multiplicity n c vs. -RoA. The solid 
line and dashed line corresponds to S/Ro — and 8/Ro = 0.2 
respectively. 



Including the multipion BE correlations, the single par- 
ticle momentum distribution for events with fixed multi- 
plicity is given in Fig. 4. It is shown that, as the multiplic- 
ity of pion increases, the number of pions at low momen- 
tum becomes larger. This has been observed in many ex- 
periments fl33| . Many different explanations of the exper- 
imental results have been proposed |34|,|35| . Some of those 
p5[ stated that the enhancement of single pion spectrum 
at low momentum spectrum is due to the two-pion Bosc- 
Einstein correlations. Our results show that the enhance- 
ment of single pion spectrum at low momentum caused 
by mult-pion correlations is more pronounced than that 
caused by two-particle correlations. So a detail analy- 
sis of the reason for the enhancement of single particle 
spectrum at low momentum is still necessary. The effects 
of the wavepacket on the pion spectrum distribution are 
also shown in the Fig. 4.. It is clear that the effects of the 
wavepacket on the pion spectrum distribution becomes 
larger for larger pion multiplicity. The single particle 
momentum distributions for total events are shown in 
Fig. 5. Similar features as figure 4 are obtained. 
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P (GeV/c) 

FIG. 4. The single particle momentum distribution for dif- 
ferent multiplicities. The wider solid line and wider dot- 
ted line corresponds to multiplicity M — 20, 5 — Ofm and 
M = 20 ,8 = 0.5/m respectively. The wider dashed line and 
wider dot-dashed line corresponds to multiplicity M — 80, 
5 = Ofm and M = 80, 5 = 0.5 fm respectively. The thin solid 
line and the thin dashed line corresponds to the input momen- 
tum distribution with 8 — Ofm and 8 = 0.5/m respectively. 
The input value of Ro and A are 5/m and 0.36GeV. 

The effects of multi-pion correlations on the two-pion 
correlation function are shown in Fig. 6, where we fixed 
the multiplicity of the events. It can clearly be seen that 
as the multiplicity of the event increases, the two-pion 
correlation function has a lower chaoticity, though the 
actual source is totally chaotic. The apparent radius de- 
rived from two-pion interferometry also becomes smaller. 
Two-pion correlation functions for different mean multi- 
plicity are shown in Fig. 7 . A similar property to that 
in Fig. 6 can be seen. It is clear that the effects of the 
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wavepacket on two-pion correlation becomes larger for 
larger multiplicity. 




P(GeV/c) 

FIG. 5. The single particle momentum distribution for dif- 
ferent no- The wider solid line and wider dotted line corre- 
sponds to no = 8, 8 — 0/m and no = 8 ,8 = 0.5/m re- 
spectively. The wider dashed line and wider dot-dashed line 
corresponds to no = 50, 8 = 0/m and no = 50, 8 — 0.5/m 
respectively. The thin solid line and thin dashed line corre- 
sponds to the input momentum distribution with 8 = 0/m 
and S — 0.5/m The input value of Ro and A are 5/m and 
0.36GeV. 




°- 8 0.05 0.1 0.15 0.2 0.25 

q( GeV/c) 

FIG. 6. Two-pion correlation function for fixed multiplic- 
ity events. The thin solid line correspond to multiplicity 
M = 2 and 8 = 0/m. The thin dashed line and thin dot- 
ted line corresponds to multiplicity M = 20,5 = 0/m and 
M = 20,5 = 0.5/m respectively. The wider solid line and 
the wider dashed line corresponds to M — 80, 8 = 0/m and 
M = 80, S = 0.5/m respectively. The input value of Ro and 
A are 3/m and 0.36GeV. 




0.05 0.1 0.15 0.2 

q (GeV/c) 

FIG. 7. Two-pion correlation function for all events. The 
solid line corresponds to no = 10, <5 = 0/m. The dashed line 
and the dotted line corresponds to no = 98,8 = 0/m and 
no = 98, 8 — 0.5 fm respectively. The input value of Ro and 
A are 5fm and 0.36GeV . 

VI. CONCLUSIONS 

Early in the next decade two heavy- ion accelerators, 
the Relativity Heavy Ion Collider(RHIC) and the Large 
Hadron Collidcr(LHC), will provide the environment for 
creating hundreds of pions. So the multi-pion Bosc- 
Einstein correlations effect in those process should be 
studied. In this paper, the multi-pion Bose-Einstein cor- 
relations and wavepacket effects on pion multiplicity dis- 
tribution, pion spectrum distribution and two-pion in- 
terferometry have been studied. It has been shown that 
multi-pion Bose-Einstein correlations and wavepacket 
cause an abundance of pions at the low momentum, in- 
crease the pion's mean multiplicity and decrease both the 
apparent radius of the source and the coherent source pa- 
rameter derived from two-pion intcrferometry. For larger 
pion multiplicity, the effects of the wavepacket on the 
pion multiplicity distribution, pion spectrum distribution 
and two-pion intcrferometry become larger. 
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The derivation of Eq. (11) will be given in this appendix. 
From Eq. (10) we have 



C 2 (pi,p 2 ) = 



{< 2|c+(p 1 )c+(p 2 )c(pb)c(pi)|2 >} 
{< l\c+{pi)c{Pi)\l >}{< l\c+(p 2 )c(p 2 )\l >}' 

(60) 



with 



\n >- 



e i P x c + 



(p)dpd 4 x) r 



10 > 



Using the relationship 

c(p)\n >= i J d i xj(x) exp(ip ■ x)\n — 1 >, 
we have 

c{p\)c{p 2 )\2 > d 4 xj(x) exp(ip ■ x) 

i J d 4 yj(y) exp(ip ■ y)\0 > ■ 
Then the two pion spectrum distribution reads 



(61) 



(62) 



(63) 



{< 2|e+(p 1 )c+(j? 2 )c(p 2 )e(p 1 )|2 >} 
{/ d^xfixx) exp(ip ■ x{) J d 4 yj*(y 1 ) exp(ip ■ y x ) 
J d 4 xj(x 2 ) exp(-ip ■ x 2 ) J d 4 yj(y 2 ) exp(-ip ■ y 2 )}. (64) 

Here {■ • ■} means the phase average(Eq.(5)). Using 
Eq.(3), Eq.(5) and Eq.(6), we have the following rela- 
tionship 

■ ■ ■ f(x n )j(yi) ■ ■ ■ j(y m )} = n^m, (65) 

then 

{j*{xi)j*{yi)j{x 2 )i(y 2 )} = {f(xi)j(x 2 )}{j*( yi )j(y 2 )} + 

{f(xi)j(V2)}{f{yt)j(x2)}. (66} 



[3] 

[4] 
[5] 
[6 

[-: 

p: 

[9 

[10 
[11 

[12 

[13 

[14 
[15 

[16 

[17 
[18 
[19 
[20 



Eq.(64) can be re-expressed as 

{< 2\c+{p 1 )c+(p 2 )c{p 2 )c{p 1 )\2 >} 
{< l|c+(ft)c(ft)|l >}{< l\c+(p 2 )c(p 2 )\l >} + 
{< l\c+(p x )c(p 2 )\l >}{< l|c+(ft)c(ft)|l >}. (67) 

Bringing Eq. (67) into Eq. (60) we obtain Eq. (11). 
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Figure Captions 

1. Pion multiplicity distribution with and without 
(solid cure) Bose-Einstcin correlations. The in- 
put value of A and no are 0.25GeV and 13. The 
dashed, dotted and dot-dashed line corresponds to 
Ro = 3/to ,8 = 0/ra; Ro = 3/m ,8 = 0.5/ra and 
Ro = 2.8/m ,8 — Ofm respectively. 

2. The mean multiplicity < M > vs. no. The solid 
line corresponds to the case without BE correla- 
tion (no =< M >). The dashed and dotted lines 
corresponds to 8 = Ofm and S = 0.5/ra respec- 
tively. The input value of Ro and A are 5/ra and 
0.36GeV. 

3. The critical multiplicity n c vs. i?oA. The solid 
line and dashed line corresponds to 8/Rq = and 
S/Ro = 0.2 respectively. 

4. The single particle momentum distribution for dif- 
ferent multiplicities. The wider solid line and wider 
dotted line corresponds to multiplicity M = 20, 
8 = Ofm and M — 20 ,8 — 0.5/ra respectively. 
The wider dashed line and wider dot-dashed line 
corresponds to multiplicity M = 80, 5 = Ofm and 
M = 80, 8 = 0.5/m respectively. The thin solid 
line and the thin dashed line corresponds to the 
input momentum distribution with 8 = Ofm and 
S = 0.5/m respectively. The input value of Rq and 
A are 5/m and 0.36GeV. 

5. The single particle momentum distribution for dif- 
ferent no- The wider solid line and wider dotted 
line corresponds to n = 8, S = Ofm and n = 8 
,5 = 0.5/m respectively. The wider dashed line 
and wider dot-dashed line corresponds to no — 50, 
S = Ofm and n = 50, 8 = 0.5 fm respectively. The 
thin solid line and thin dashed line corresponds to 
the input momentum distribution with 8 = Ofm 
and 8 = 0.5 fm The input value of Ro and A are 
5/m and 0.36Gey. 

6. Two-pion correlation function for fixed multiplicity 
events. The thin solid line correspond to multi- 
plicity M = 2 and 8 = Ofm. The thin dashed 
line and thin dotted line corresponds to multiplic- 
ity M = 20,(5 = 0/m and M = 20,8 = 0.5/m 
respectively. The wider solid line and the wider 
dashed line corresponds to M = 80, 8 = Ofm and 
M = 80, 8 = 0.5/m respectively. The input value 
of Ro and A are 3/m and OMGeV. 

7. Two-pion correlation function for all events. The 
solid line corresponds to no = 10, 8 — Ofm. The 
dashed line and the dotted line corresponds to no — 
98, 8 = Ofm and n = 98, 8 = 0.5 fm respectively. 
The input value of Ro and A are 5/m and 0.36GeV. 
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